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Abstract 

The effect of quasi-particle (QP) 'scattering' by the vortex lattice on the 
de-Haas van-Alphen oscihations in a pure type-II superconductor is investi- 
gated within mean field, asymptotic perturbation theory. Using a 2D electron 
gas model it is shown that, due to a strict phase coherence in the many- 
particle correlation functions, the 'scattering' effect in the asymptotic limit 
Ep /hujc ^ 1) is much weaker than what is predicted by the random vortex 
lattice model proposed by Maki and Stephen, which destroys this coherence . 
The coherent many particle configuration is a collinear array of many particle 
coordinates , localized within a spatial region with size of the order of the 
magnetic length. The amplitude of the magnetization oscillations is sharply 
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damped just below Hc2 because of strong 180° out of phase magnetic oscilla- 
tions in the superconducting condensation energy , which tend to cancel the 
normal electron oscillations. Within the ideal 2D model used it is found, how- 
ever, that because of the relative smallness of the quartic and higher order 
terms in the expansion , the oscillations amplitude at lower fields does not 
really damp to zero, but only reverses sign and remains virtually undamped 
well below Hc2- This conclusion may be changed if disorder in the vortex 
lattice, or vortex lines motion will be taken into account. The reduced QP 
'scattering' effect may be responsible for the apparent crossover from a strong 
damping of the dHvA oscillations just below Hc2 to a weaker damping at 
lower fields observed experimentally in several 3D superconductors. 



PACS numbers: 74.60.-w, 71.25.Hc 

I. INTRODUCTION 

Magnetic Quantum Oscillations have been recently observed in several type-II supercon- 
ductors below Hc2 IH]- [§. A systematic study of this remarkable effect has been impaired, 
however, by the lack of a complete quantitative theory of the de Haas-van Alphen (dHvA) 
effect in the vortex state, analogous to the Lifshitz-Kosevich (LK) theory in normal metals 
0. Such a theory would require a detailed analysis of the effect of the superconducting 
order parameter on the magnetization oscillations in the vortex state, which turns out to be 
an extremely subtle theoretical problem. 

A common feature reported by all experimental groups so far,which is far from being 
well understood, has been the observation of an additional damping in the dHvA amplitude 
below Hc2- Several theoretical papers have attributed this attenuation to the broadening 
of the Landau levels by the inhomogeneous pair potential, building up below Hc2- There 
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is, however, a remarkable disagreement among the various theoretical approaches to this 
problem concerning both the size of the attenuation factor and its detailed dependence on 
the strength of the field below Hc2- 



The semiclassical approach, adopted originally by Maki , elaborated later by Stephen 
TT| , and reviewed very recently by Wasserman and Springford []T^ , considered the correction 
to the Quasi-Particles (QP) lifetime due to the 'scattering' by the vortex lattice. 

For the first harmonic of the oscillatory part of the magnetization below Hc2 they pre- 
dicted: 

Mosc = Mn,osc exp [-A (Ao, np)] (1) 

with 



A(Ao,n^) = vr=^/^^ (2) 



rip 

where Mn,osc is the corresponding normal electrons contribution to the oscillatory magne- 
tization, Ao = Ao/hujc, Aq - the magnitude of the superconducting (SC) order parameter, 
Uc-the cyclotron frequency, and rip = Ep/Tiojc is the Landau level index corresponding to 
the extremal orbit on the Fermi surface. 

An exponential damping with a different exponent, i.e. A ~ Ao/n]/^, has been proposed 
by Norman et al. [|l^], who carried out a full quantum mechanical calculation, based on a 
numerical solution of the Bogoliubov-de Gennes (BdG) equations for the quasi particles at 
low temperatures. The numerical computations carried out by these authors were limited, 
however, to relatively small values of ^/np- 



A similar approach, invoked by Dukan and Tesanovic ||T^, has led to a qualitatively 
different behavior at very low temperature, that is a power law attenuation of the dHvA 
oscillations below Hc2 

Mosc/Mn,osc ~ {kBT/ Aof + [{kBT/ Ao)'] (3) 

This result was obtained by considering only the dominant contribution to the dHvA 
oscillations as originating in the 'gapless' region of the QP spectrum around the Fermi 



surface. It looks like a low temperature high A expansion, emphasizing the opening of the 
SC gap well below ifc2 

In the Maki-Stephen (MS) theory it is assumed that the vortex lattice acts like a random 
potential for the quasi particles, and so by averaging over the realizations of the vortex lattice, 
the QP self energy acquires a large imaginary part, leading to a strong exponential damping 
of the dHvA amplitude. 

Strictly speaking, however, the 'scattering' by the inhomogeneous pair potential is a 
highly coherent process, as in multiple Andreev reflection at the interfaces of a 2D periodic 
array of normal and SC phases [ p^ . 

Thus, in the ideal, self consistent vortex lattice model the QP self energy has zero 
imaginary part. However, the broadening of the Landau levels into real energy bands by the 
inhomogeneous pair potential should lead to damping of the dHvA oscillations even for an 
infinite QP 'lifetime'. The term 'scattering' refers to this inhomogeneous broadening effect 
in the present paper. 

The quadratic dependence of X^AqjUf) on Aq in Eq.( 0) reflects its origin in a pertur- 
bation expansion of the the QP self energy in the mean square order parameter, which is 
strictly valid for a pure type-II superconductor only at sufficiently high temperatures, when 
the Lifshitz-Kosevich § thermal smearing factor X = 271"^ kBT/huc > 1. 

It is therefore very interesting to compare the result expressed in Eq.(|l]) to that obtained 



by Maniv et al. |lT6| , p!7| , who considered Gorkov's expansion of the SC free energy in the 
small vortex state order parameter near Hc2- 

Using a semi classical approximation, valid for y^np ^ 1, these authors have found a 
quadratic term consistent with Eq.(|]), but a quartic term M^f^/M„, 

osc ~ ^o/'^F ) which 

is smaller than the corresponding term, obtained in MS theory (i.e. Mj^jl/ Mn^osc ~ ~ 
Aq/uf), by the factor l/^yrTp <^ 1. 

Very recently, Bruun et al. |T8| developed an exact numerical scheme for calculating the 
coefficients of the Gorkov's expansion within the same model used by Maniv et al. and 
found good agreement with the results obtained by Norman et al. for small values of np. 



They have also made an estimate of the np dependence of the quartic term by using an 
approximation similar in spirit to the random lattice approximation, and found a result 
which agrees with that of Maki and Stephen. 

In the light of this controversy our purpose in the present paper is to carefully examine 
the high temperature X > 1, small Aq, asymptotic {^/np ^ 1) expansion, in order to 
elucidate the origin of the disputed np dependence. We find that, incoherent 'scattering' 
channels, which generate the dominant contribution in the random lattice approximation, are 
completely cancelled in the self-consistent , periodic lattice calculation, due to the presence 
of a strict phase coherence in the four-particle correlation function. The remaining coherent 
four-particle configuration is a coUinear array of four-electron coordinates, localized within 
a spatial region with size of the order of the magnetic length. 

As a result, the inhomogeneous broadening of the Landau levels by the pair-potential 
does not contribute significantly to the damping of the dHvA amplitude just below Hc2 
, as is the case in the MS theory. The dominant damping mechanism in the asymptotic 
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limit ,np ^ 1 is found to arise from the strong, 180° out of phase oscillations of the SC 
condensation energy with respect to the normal electrons oscillations, as was first proposed 
in Ref . . 

The organization of the paper is as follows: In Sec. II we present the general framework of 
our approach, which is based on the Gorkov's expansion of the SC free energy near Hc2- In 
Sec. Ill we carefully examine the quartic term in the asymptotic limit i/rTp ^ 1 and justify 
the main approximation used in our calculation. In Sec. IV we derive simple analytical 
expressions for the self consistent order parameter and for the oscillatory magnetization , 
and verify the validity of our claimed np dependence in the asymptotic limit. In Sec.V we 
discuss the connection of our theory to the other theoretical approaches and compare our 
predictions to experiment. 
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II. SMALL ORDER PARAMETER EXPANSION 



We consider the quadratic F^^^ and the quartic Fj^^ terms in the Gorkov-Ginsburg- 
Landau expansion of the SC free energy in the SC order parameter A (1^) For the 

sake of simphcity let us assume a 2D electron gas model and neglect the spin degrees of 
freedom. The former assumption may be justified by noting that the main contribution to 
the dHvA effect in an isotropic 3D normal electron system comes from the extremal orbit, 
corresponding to the value kz = of the electron momentum parallel to the field direction, 
and to the Landau level n ^ np = 

Note, however, that the pairing effect responsible for the Cooper instability in a 3D 
electron gas is dominated by the region near the Fermi 

surface. Since the focus in the present paper is on the effect of QP 'scattering' by the 
vortex lattice, we shall ignore, for the sake of simplicity, this aspect of 3D systems in our 
calculation, but will return to this problem later while discussing our results in connection 
with experiment. 

Thus 

Ff) = ^Jd^r I A (r)l' - J £nd\2K2 (^1, 7^2) A (^i) A^ (^2) (4) 

and 

Fi^) = IJ d'nd\2d'rsd\,K, ({r,}) A (r^i) A^ (^2) A (r^s) A* (^4) (5) 

where V is BCS interaction constant. The kernels K2 and K4 are expressed through 
the product of two and four electron Green's functions Gq (^'f^ 1,^1+1, {—Y ^uj, / = 1,2 
with "P's = "P'l, and / = 1,...,4 with T's = Y^i, respectively, in magnetic field (in 
our case^ = (0,0, iJ)). Here cj^, is the Matzubara frequency. The magnetic field 
breaks translational symmetry of BCS Hamiltonian and Green's function. However, ow- 
ing to the gauge symmetry the Green's function can be represented in a factorized form: 



Go {T^i, T^i+ijUJiy) = g {Y^i, 'r^i+i) Gq {puuJu) |T9[, where the reduced Green's function 
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depends only on relative coordinates ~ffi = {'r'l+i — 'r'l) /an, = 0, o,h = \f^- 

this expression L„ is the Lagaurre polynomial of order n, and /i is the chemical poten- 
tial. The non invariant gauge factor in the symmetric gauge has the form: g {l^i, 'r^i+i) = 
exp {-iEik {ri^i + n+i^j) pi^k/^a-n), £ik = Ski is antisymmetric tensor in 2D space. 

In quasiclassical limit ( ^/r^F ^ 1 ) Green's function Gq {pi, u,) has two different types of 
behavior. Near turning point pi c:^ 2rF { rp = \p2jip is the cyclotron radius ) it is a smooth 
function, and in the intermediate region pi <2rp it is a sharply oscillating function. 

From now on we shall express all spatial variables, except for l^j, in units of the magnetic 
length. 

Using asymptotic for Lagaurre polynomial L„ at n ^ oo pO| and the Poisson summation 
formula for the sum over n in one can show that Go (p;,u;^) ~ V'^-f^ 4r|, — pf < 
{IQnp 1^)^^^ , and 



ie,.,..J (uj,A exp 



iSu^, [np (0 + sin0)] - 0^ 



(^o{pi,UJu) - -——rw—r, TJT—r, ^^74 V) 



where J ^ (tu^) = 1 - exp {2'Rie^^np - 2ix \ujy\ juj^), sin0 = f-(\- [^Y^ 
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sign (uj^), at pi < 2rp. Note that this expression is a generalization of the Green's function 



obtained in [16| for pi/2rp <^ 1. The similarity to the short distance limit is apparent since 
the analytic continuation of the Lagaurre polynomial asymptotic in the complex n plane 
gives rise to short distance approximation at \n\ — > oo. 

Combining the gauge factors we obtain the dependence of K2 and on the two and 
four-particle center of mass coordinates = {^^1=1^1^ /'^o-h and R = /4:aH 
respectively: 

K2 (T^i, y'a) = exp {-ieikTipk) ^ K2,u (p) (8) 
^4 ({^i}) = exp {-^leikRiQk) - E ^4,. ({pj) (9) 

P V 



Where = = - {"^i - + - ^ a) /§ and 

K2,u (p) = Go (p, -Uu) Go (p, CO,) (10) 

-^4,1. {{Pl}) = Go (pi, -UJ,) Go (P2, t^,.) <5o (P3, -t^i.) Go (P4, t^i/) (11) 

For the order parameter we use the Abrikosov form. In symmetric gauge 



A ("r*) = Ao exp (— Z/^ + ixy^ ^ exp 



, , 7™\^ nn , 



:i2) 



The numbers 7 and a^; are arbitrary. It can be shown that this quasiperiodical form results 
from gauge symmetry as a solution minimizing free energy [|l^. The order parameter is 
normalized as 



where is the number of vortices. 

Using expression (|1^) for the order parameter in Eq.(^ for the quadratic term we get 
after integrating over the two-particle center of mass coordinate 'r* 

Fi^) ={\^-a) Tra^NAl (14) 

where A oc 4 E / ^^^£""'^^^2,^ (p)- The gaussian factor restricts the effective integration 
region by the distances p ~ 1. This important fact means the loss of coherence in electron 
pair propagation over distances much larger than the magnetic length. 

The localization of the electron correlation functions within a region of a size of the order 
of the order of the magnetic length can be shown to exist, in quasiclassical limit, also for 
many electron configurations. Here we discuss only four particle correlations. Integration of 
(1) over R gives rise to expression 

Fi^) = / d'Qe-'<^' J d'Sd'TV ( S^, 7^) i ^ K.,. {{pi}) (15) 

where = ^ (7^3 — T^i), ^ = i (7^4 ~ 7^2)- The function V , 7^^ includes the summa- 
tion over vortices: 



mi ,m,2 



exp 



- (mi + 2SyY - (m2 + 2Tyf (16) 



For convenience, we denoted mi = it {n2 — ni) /a^, m2 = vr (723 — rii) /a^, where Ui is the 
summation index corresponding to the order parameter A (l^i) in Eg. (p!2D . Note that only 
a 2D sum remains in Eq.(|16D from the original 4D sum; the free double sum is equal to the 
number of vortices- A^, while the summation in Eq.(|l6D is invariant under the vortex number 
shift transformation Ui ^ rii + uq. Note also that, the combination of gauge factors from 
the order parameters and Green's function causes integral over center mass coordinate R 
to vanish if rii — n2 + — ^ 0. 

It is clear from Eq.(^B[) that the lattice sum in V , T^^ is dominated by the lattice 
point mi ~ —'2Sy, m2 — —2Ty. Taking only this term in the sum into account we obtain 
that 

V ( S^, T^) oc A^e~''*^^'^'^^'^«+^^'^"+^"'^''^ (17) 



The lattice function V ,T^^ being multiplied by the kernel K4^i_,{{pi}) determines the 
free energy distribution in the space of relative electron coordinates. 

III. ASYMPTOTIC LOCAL APPROXIMATION 

In this section we analyze in detail the relative importance of all different spatial regions 
contributing to the multiple integral , Eq.(|r^) , which determines the quartic term in the 
asymptotic limit ^/nF ^ 1. 

Every four electron configuration in this integral is defined by the three vectors 
(^Q,^,^^. Since the kernel K4^,^(^Q,^,^^ is a bounded function at Q » 1, we 
may conclude from Eg . (pTsD that the main contribution to the free energy comes from the 
the region Q < 1 , and so assume in what follows that Q < 1. To estimate the integrals 
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over and we separate the entire domain of integrations into three regions according to 
the behavior of the Green's function (Eq.(^)) and the function V , T^^. 

1. The turning point region pi ~ 2rp, for which in the asymptotic hmit y^np ^ 1, 



Pl,3 



P2,4 



5 ± 5^ -2S -2rF = 2v/2^ > 1 



2 1^ ^ 7^ ~ 2T ~ 2ri. = 2^/2^ > 1, 



(18) 



The size of this region is of the order Api = AS = AT ~ 1/r]/^. It is characterized 

~^ — 1/3 

by a smooth behavior of all Green's functions Go (Pi) ~ GQ{2rp) ~ 1/''^^? and kernel 



i^4,. ~ [Go (2rp)] r^l/nf. 

2. In the intermediate region 1 <^ pi < 2rp variables S and T are still large: S ^ 1, 
T ^ 1 but essentially less then rp. Here the Green's function is approximated by the 
expression Eq.(0). To simplify the considerations we note that the oscillating phase factor 
in Eq.(^) can be written as 



Hp ((f) + sin( 



(pf^rp) p 



(19) 



where n/A < ^ (x) < 1 at < a; < 1. For our purpose we may take ^ {x) = 1. Substituting 
Pi 3 ~ 2S, p2,4 — 2T, the kernel ^4 in this region can be transformed to 

i4^/2E^{S-T) 



exp 



oc 



ST [rl-S^f^ {rl-T^) 



1/2 



(20) 



3. In the third region all variables Q, are of the order of the magnetic length: 

Q,S,T ^ 1. Here P(5^,f^) becomes a smooth function and the strong oscillations of the 
integrand in Eq.(|T5|) arise only from K^^^. Since ^ (x) — 1 at x — the exact short distance 
asymptotic is 

1 



Ki,u ({Pi}) oc 



exp 



zV2rip (pi - P2 + P3 - Pi] 



(21) 



np (P1P2P3P4)^^^ 

Let us first consider the case when both variables belong to the first region (i.e. 

region(l,l)): Replacing K/^^y (^Q , T^^ by a constant of the order ~ ^/n^p^ we can estimate 
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the free energy ,Eq.(^), for a given Matsubara frequency in polar coordinate system 5^ 
{S,as), = (T, ctf) as 



rp 2iT 



Ptu^-JJ^ J SdSTdr j dasdate-^'^^'*'^''^^'''^ (22) 



'^F rp-Ap 



where (a^, at) = 7 cos (a^ — at) + sin (a^ + at) — 7 cos (a^ + at). After simple integration 
over as and at by the stationary phase method, which is justified by the very large values 
of ST, Eq.(^) transforms to 

1 ''^ 1 
F(fJ^4j^ f dSdTe-'^^^^^ ^ (23) 

where (f)f = a/1 + 7^ ± 7 are the values of (j) (os, at the stationary points. 

Let us next consider the case when both S and T belong to the intermediate region (i.e. 
region (2,2)): Here K^^^dpi}), defined by Eq.d^) like ©(S^, 7^), is a sharply oscillating 



function with frequency of order of y/np. Substituting (0),(|r3) into (|1^) we obtain for 
1< 5,T < 

d^Sd^Texp [4i (^27^ (S-T)- 2-fSyTy - SyT^ - S^Ty 



i?{4) 



OC 



(24) 



ST{rl-S^f'^ (4-T2)'/' 
The integrals over angle variables are identical to previously considered ones and give a factor 
of the order 1/ST. Thus free energy (pi) reduces to F^^} oc f '^^'^'^'^^^I'^^i^^^^j-^ '^^='''^^^0] _ 
Taking for the smooth preexponential factor its value at a point S ~ ^/np.| T ~ ^/np and 
noting that there is no stationary point for the integrals over S and T, we get the following 
estimate for free energy in the intermediate region: F^'^^ ~ 1/n^p. 

It is clear that if one of variables, for example S", is from the turning point region and 
the other, T, from the intermediate region (i.e. region (1,2)), the free energy F'^'^^ will be 
proportional to a factor of the order ~ l/v?p^. This result follows from the fact that the 
product of two of the Green's function are proportional to ~ l/rip'^ , the product of the other 
two is proportional to ~ l/np and the integration over S and T yields the factor ~ l/np. 
The integration over the angles as and at produces the factor ~ 1/ST which cancels the 
same factor in the numerator of the integrand (pTSf). 
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Let us consider now the four-electron configuration when one of the variables, S or T, 
e.g. T , is of the order of the magnetic length T ~ 1 (i.e. region (3)). In contrast to previous 
cases, S,T ^ 1, where the oscillations with Q can be neglected, in this case an additional 
small factor arises from the integration over Q . It is partially cancelled by the large factor 
arising from the integration over ^ if it is from the turning point region. Assuming that 
S ^ rp and T ~ 1 (i.e. region (1,3)), the estimate of the free energy (|15|) is given by the 
formula 



exp 



(25) 

(26) 



The integral over Q can be divided into two regions: Q < T and Q > T. The main 
contribution comes from Q < T, where the above exponential factor does not depend on the 
component Q", parallel to . For the small transverse component we get 



+ + f ~2T+(- + -) (q^)'' 

\P2 pi) ^ ^ 



(27) 



Thus the integral over is proportional to ~ I/'^-f^ ^^"^ ^^e integral over gives (iQ" ~ 
T. Performing the integration over angles and at gives, as previously, 1/ST and taking 
into account that integral over S from a smooth function can be estimated by the area of 

1/3 

the integration region: SAS ~ rpAp ~ Up , we get for Q <T that 

1 

n'p' " •> " " ' ' n 



^iJ ~ J TdTexp {-At^T (l ± 



31/12 



(28) 



In the second region ,Q > T, the integral over Q leads to a contribution smaller by the 
factor n]/^ than (^). This is because the phase (p7|) does not depend on and we can 
neglect the quadratic in T-*- term in comparison with the linear one in Eq.(^).Thus the 
integral over Q can be estimated here as ~ 1/np. 

If belongs to the intermediate region (i.e. considering the region (2,3)), the contribu- 
tion to the free energy will be even smaller because of the Green's function oscillations and 
smaller preexponential factor. 
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Finally, in the short distance region Pi < 1 the smooth function V , 7^^ is of the order 



one and we should calculate the integrals in (15) from the kernel {{pi}) Eg. (pT]) . The 
phase factor of K^^^ which is proportional to 

$ S^, T^) = pi - p2 + P3 - P4 = 

= 2{|^+ S^l-I'^- T^l + I'^- 5^|-|^+ T^|} (29) 

at fixed Q has a set of stationary points S < Q, T < Q ioi collinear vectors Q , Going 
back to variables we conclude that this configuration is equivalent to the propagation of 
"odd" particles in a single direction ft and "even" particles in the opposite direction: 



P 



Pin, p 2 



-p2n, P3 = psTi, p 



-pin, 



(30) 



where it is an arbitrary unit vector. Note that the special role of the configuration ( pOD 
follows from the fact that for this configuration ']} ^ = $ (pj) 7f = in the integration 
region. 

It should be emphasized that in the resulting coherent configuration the correlation 
among all four particles is essential. The phase factor, Eg . (|29|) , and stationary point equa- 
tions can not be factorized. 

Expanding ^ {^Q , S^, in the coordinates S-^ and T-*- perpendicular to ft , and noting 
that from the definition of Q follows it = —Q /Q, we can reduce (|29| ) to 



.Pi ^ P3, 



where Pi + P3 = P2 + Pi = ^Q- Now performing the integration of the kernel 7^4 ^, ({pi}) , 
given in Eq.(^), over S-^ and T-*- we obtain a factor ~ (pip2P3P4) ' ^ which gives rise to 
the V dependent free energy ~ l/n^p"^. 

This result completes our analysis: it implies that in the asymptotic limit ,^Jrip ^ 1, the 
dominant contributions to the quartic term of the free energy originate in the short distance 
region only. Thus, using the short distance approximation for Ki^^y, taking into account its 
dependence on ujy, and restoring the exact form of T> , T^^ ,Eq.(^^, the quartic term of 
the SC free energy can be written in the form: 
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exp [-4(52 + T2)sin2^ + 4z5Tsin(2^)] ^ exp 



(32) 



mi ,m2 



—2i'~fmi'm2 + 4z {rriiT + 77225) cos — (miS* + m2T) sin 0] 

where the angle 6 describes the direction of the unit vector it , 2qi, = J {—ojy) J (u^), and 
= 2(2z/ + l)aH/(, with ^ = Tivp I ij^kBT). 

This free energy is determined by colhnear, essentially four particles configurations with 
the size of the order of magnetic length. The resulting expression, Eg. (P^) , agrees with the 
quartic term Fj^^ derived previously by Maniv et al. |rB|. Our present considerations justify 
the used approximation and clarify the geometry of coherent configurations. 

It should be emphasized that the random vortex lattice approximation, used in , gives 
rise to a markedly different result, namely F^^} ~ \jnp. The reason for the disagreement 



is due to the averaging over random vortex lattice |TT[ which leads to factorization of the 
multiple products of pair potentials into products of pair correlation functions only. 
For example, the quartic term in this approximation becomes 

(A (r^i) (^2) A (r^a) A^ (^4)) oc (A (r^i) A^ (^2)) (A (r^s) A^ (^4)) 



oc exp 



exp H (C (2,1) + C(4,3))] 



where C (2, 1) = (x2 + X\) (7/2 — Vx) is the (Landau) gauge factor of G'o(^i5 "^2; '-^v)-, and 
(....) stands for averaging over vortex distributions. The corresponding free energy is given 
by 

1 



Pf) « / n^'^*^4,.({p^})exp 



2 I 2 
Pi +P3 



X 



exp 



- [(0:4 - X2) {ys - yi) - (a;3 - xi) {y^ - 2/2)] 



(33) 



In the important region of integration pi,p3 < 1, the constraint J2i=i l^i = 0, implies that 
7?2 ~ ~"P^4) oi' alternatively T^i ~ 1^2 and ~ T^4. Thus the imaginary exponent in 
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Eq.(p3D is always of the order unity or smaller in the important region of integrations, and 
so the total gauge factor is a smooth function. Consequently, only two distances in the 
kernel Ki^^y are restricted to the size of the order of magnetic length, allowing two others 
to be arbitrary. Thus, in contrast to our result, where all long range configurations of 
electron pairs in the ordered vortex lattice interfere destructively (i.e. appear as incoherent 
'scattering' channels) , the smoothing of rapidly oscillating gauge factors in the MS theory 
introduces a huge 'incoherent' contribution to the SC free energy. 

Let us estimate now the quartic term within this approximation. Substituting for the 
Green's functions Go (pi, —ijJv) and Go (ps, —^u) their approximants in the short distance 
region, and omiting the smooth gauge factor, we get for the free energy -Fj^'', after integrating 
over the center mass coordinates: 



1 



oc 



n 

exp 



Tj2^ \Y.y l] '^^^ (/^2, UJy) Go (p4, UJu) 



1/2 / / 

F [pips) \l=l 

(i\/2n^(pi + P3jj exp 



(34) 



The main contribution to Eq.(|3^ arises from the region p2 ^ 2rp, ^ 2rp. Allowing 
7?2! 7^4 to vary independently within the turning point region, "p^^, 7^3 are not independent 
variables; taking 7?^ as the third independent variable of integration , and noting that 
7?i — —7^3, we have 



1/2 
Up 



(fpl- 



Pi 



(35) 



Now, since / GQ{p)d p ~ 1, and the integration over Jf^ yields the factor ~ l/i/Snp, 



the resulting np dependence is 1/np, in agreement with [p 



IV. SELF CONSISTENT ORDER PARAMETER 

The local approximation, verified in the previous section, becomes very transparent if we 
rewrite the free energy ,Eq.(^), as a functional of the order parameter A (7^) . After some 
straightforward, but combersome calculations one can show that Eg . (|3^) is equivalent to 
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) oc^ J(fRj dOdQf (Q) e-^Q' / / ( + (5 + T) Tf ) (36) 
A* + (5 - T) rf ) A - (S + T) It) A* ['^-{S-T)- 



where f (Q) = Eg^e-^"'^^, and A (^) is defined by Eq.(|12D. Since \S\ , |T| < Q < 1, the 
expression (|36| ) can be considered as averaging of the four order parameter product over a 
region with radius of the order of the magnetic length. The additional averaging over the 
direction of rf in (|32|) leads to a completely local expression plus a nonlocal correction, i.e.: 

F^'^ = Bjd^R\A{-n)\' + F^X (37) 

where B oc -^72 I dQf (Q) e'^'^' J dS J dTe~^^^'^^'\ The nonlocal correction , F^tL, 
is numerically small since it arises from high (i.e fourth and higher) order terms in the 
'cumulant' expansion of the exponential in Eq.(p2D (see Ref. p^). 

This result is of fundamental importance since it shows that the well known, fully local 
form of the Ginzburg-Landau free energy functional in the low field regime near Tc{H = 0) 
is basically valid also in the opposite , high magnetic field regime near Hc2{T = 0). This 
locality is closely related to the coherence effect discussed above. For example, in the random 
lattice approximation, discussed in the previous section, the dominant contribution to F^^^ 
is extremely nonlocal. 

Neglecting the small nonlocal correction, the total SC free energy, up to fourth order in 



Aq, can be turned into the following one parameter variational form 

B 



F, 

''H 

-2 



fs = -T-j- 2~ ~ -^21) 

Nnai 



(38) 



where 'D2D = wic/Svr^ (i.e. the 2D single electron density of states), 

a = 2^Y.M<luhu-l/9 (39) 



with 7^ = j;^ dpe-^'^P-^p'', g = VV2D , and vd = {Td/2T - 1), where To is the Debye 
temperature. 
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The coefficient ,B, of the quartic term can be readily obtained from Eq.(p7D (after re- 
placing Q with p/2): 

B=p^^{^rj:R<Qi)s. (40) 

with 61, = 2n /q°° dpe~'^°'''P~''^erp{p/y/2), = hvp/TrkBTc, and Pa is the geometrical factor 
of the Abrikosov lattice [|T^ . 

The key parameters , which control the crossover from the low field to the high field 
regime, are an/Cy ^ = ^n'^kBT/huJc] they are connected by: 

X = 2vr(2np)i/2(^) (41) 

which means that in the asymptotic limit considered here, our 'high' temperatue regime, 
X ~ 1, of the quantum magnetic oscillations domain is still in the low temperature regime 
of the SC- normal phase boundary, since ^ ~ l/y/np ^ 1 

In this case ^ 2 for all z/, and the coefficient, a, of the quadratic term, can be 
calculated from Eq.(pP|) by dividing the sum over the Matzubara frequencies ,z/, into two 
regions: (1) a,, <^ 1, namely u <^ i^max = ^/2"\/2, and (2) u > Vmax- The contribution 
from the ffist region is /q^ e^^(l — erf{x))dx ~ 1.147, while the sum in the second region 
leads to the famihar logarithmic expression YX=v^ax ^1^^ + -'-Z^) ~ In (v^-^^^l )^ provided 
that the Debye cut-off temperatue To = {2i'o + 1)T is much larger than {2i'max + ^)T. The 
last condition may be rewritten in a more transparent form, i.e. (kBT^i/tiuJc)'^ ^ np/2T[^. 

Combining the contributions from the two regions we find 



a f» In 



(42) 



.v^e(o). 

where ^(0) = .IShvp/kBTc ^ .56^o- 

Now consider the coefficient B, of the quartic term. Again, we divide the Matzubara 
sum into the same two regions: In the ffist, where u <^ i^max, each term,(5y, is indepen- 
dent of u so that Ylu=1f ~ 4:^max Io° dpG~^^ P^dp ~ •63(^), whereas the second region 

r 1 1^ 1 / 

yields J2'^=u,^^^ 2{2u+i)aH/(; ~ w'^/^h- Combining these results we find that the sum over 
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Matzubara frequencies changes significantly tlie np dependence of the quartic term with 



respect to the indevidual Fj^^ terms, since Y.v^v ~ -^^i 



C ^ 1/2 



an • 



Ur 



We thus find that B/inkBT^f ^ 1.38/ EpfiuJc, so that 



fs 



(43) 



It should be emphasized, here again, that the np dependence of both the quartic and the 
quadratic terms in fg above differs by the large factor ^/nF from the indevidual terms 
because of the sum over the Matzubara frequencies. 
Using Expression (^31) , the self consistent mean field order parameter is given by 

This expression is identical to the well known high field limit of the Gorkov-Ginzburg- 
Landau SC order parameter ( ||12[)- Indeed, at magnetic fields H near Hc2{0) = 0o/2vr^(O)^, 
00 = ch/2e, where an ~ v^^(O), we have: 

, so that Eq.(PD reduces to the well known result 

1/2 



Ao ^ 1.7kBT, 



inf^y 

. \V2mJ 



1.7kBT,[l~H/H,2{0)f^ (45) 



Interestingly, the np dependence of the self consistent Aq obtained in Eq. (^4D for H ^ Hc2 
determines a small parameter: 

^ = ^ ~ .36 [1 - H/H,2m 
n p 

which is seen to be the expansion parameter in the perturbation theory leading to Eq.(^). 
This observation will be further discussed in the next section. 

Note that in deriving the above expressions for the self consistent order parameter we 
have neglected the oscillatory compomemts of the SC free energy, which should add an 



oscillatory contribution to the order parameter [|T^,|l^]. In the 'high temperature' limit 



considered, this oscillatoray term is much smaller than the nonoscillatory one, except for a 
very narrow region near H^.^ ||27|| . 



Let us consider now the magnetization oscillations; the dominant contribution to the 
superconducting part can be obtained by differentiating the density of states factors in 
the free energy with respect to magnetic field, namely: 



M,„,,oc-V^^ (46) 



Explicitly we have: 



irkBic C 



For X > 1, ^ ^ —{Snnp/ H) sin {2TTnF)e~^'^'^^^^^ , so that the sum over u is limited by 
the thermal damping factor to the first few terms only. This contrasts the nonoscillatory 
magnetization, which picks up contributions from many Matzubara frequencies. 

Thus the first harmonic of the oscillatory magnetization, Mosc, just below Hc2 can be 



written as 



M = —M ^ M 



F 

X 



1/2 ' 3/2 

Up 



(48) 



where (3a ~ 1-16 for a triangular lattice , and M^^osc = -pr sin (27rnF) is the normal electrons 
oscillatory magnetization . 



In the expansion ( ^81) there are two scales of order parameter Aq. Near Hc2, where 
A2 < {hucf/^El^^ (i.e. Al < which means that In [ y/||(o) ] ~ -'-/^f^ )' attenuation 
of the magnetization oscillations amplitude occurs as the result of the electron pairing. Here 
the contribution of the many electron coherent configurations are negligeable. Far from Hc2, 
where Aq ~ TiuJcEp (i.e. ^ np so that In [-^7=|^] ~ 1 ) , the quadratic and the quartic 
terms in the free energy (and magnetization) are comparable. It can be shown that 
the higher order terms in this expansion are determined by the parameter Aq/TiUcEp = 
In the region where this parameter is of the order unity or larger the SC state is a highly 
correlated many electron-pair configuration , which is quite different from the condensate of 
electron pairs, dominating the SC free energy just below Hc2- 
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V. CONCLUSION 



The results of the last two sections enable us now to critically discuss the various the- 
oretical approaches to the problem of the intrinsic attenuation of the dHvA oscillations in 
the vortex state, and the relevance of our model to real experiemts. 

It is, first of all, clear that the assumption of disordered vortex lattice, and the conse- 
quent averaging over the random pair potential configurations, which greatly simplified the 
analysis in the MS theory [|10|, [|lT| , replaces the many electron correlation function with 



a product of pair correlation functions, and so greatly overestimates the QP 'scattering' 

1/2 ~ 

effect in the asymptotic limit Up ^ 1 . In fact, up to the second order in Aq, our result 
(Eg. (|48|)) is identical to that obtained by MS (Eq.(|l|)). The higher order terms, however, 

1 /2 

differ substantially; our quartic term is 1/np -C 1 smaller than that obtained by expanding 
the exponential in Eq.((|I])) up to second order in A(Ao,r;,i;'). 

This result reflects a very interesting phenomenon: In the ground Landau level approxi- 
mation for the condensate of Cooper-pairs, the quadratic term in the free energy expansion 



is known [|l^] to be completely independent of the vortex lines distribution. Therefore, it has 
nothing to do with the broadening of the Landau levels by the inhomogeneous pair potential 
in the vortex state. Indeed, in the standard expression for the SC free energy in terms 
of the 'dressed' electron Green's function (or the QP Green's function) the entire series of 
self energy corrections is multiplied by a second order factor in A. Consequently, the quartic 
term is the lowest order correction to the free energy, which contains the 'scattering' effect. 
It may be, therefore, concluded that in the asymptotic limit of the 2D model used here, the 
'scattering' effect is much weaker than what predicted by any theory consistent with the 



random vortex lattice approximation |]TT| |T0| ||18|| . 

The structure of our expression for the free energy , (^3|) ,as well as for the oscillatory 



magnetization, (^), suggests that the small expansion parameter in the theory is x 
rather than as suggested by Eq.(|l]). The full expansion should therefore read: 

"l - 7i^/^y/^xe{x)] (49) 
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M ^ M 



where at x <^ 1 the function (x) has an expansion 

Now the expression within the square brackets in Eq.(^) vanishes at xQ{x) = 1/ ir^^'^n'J'^ . 
Thus a sign inversion of the magnetic oscillations amplitude takes place at a; ~ 1 / rr'^^'^np'^ <^ 
1, where 0(x) ~ 1, i.e. well within the range of validity of our expansion p!3| , |23| . 

One therefore expects that in a 2D superconductor the dHvA amplitude will reverse sign 
due to pairing at a certain field ,i/j„i,, below Hc2, and remains virtually undamped well 
below the point of inversion. 

This conclusion may be changed if disorder in the vortex lattice, or vortex lines motion 
will be taken into account, as indicated by the MS result. However, application of the MS 
model to real disordered vortex lattices should be considered very cautiously since the effect 
of disorder has not been introduced self consistently there. 

The crossover to the low temperature power law behavior, obtained by Dukan et al. , is 
reflected in our theory by the breakdown of perturbation theory at very low temperature. 
At such low temperatures, the LK thermal smearing parameter X <^ 1, and our expansion 

does not exist for all magnetic fields since the density of states parameter q^y diverges like: 

2 2 
" {l-exp[(2z/+l)X]} ^ (2z/+ 1)^X2 

when a Landau level crosses the Fermi energy with half integer filling factor np. 

Under this condition, and for sufficiently small Aq, the SC pairing is restricted to a single 
Landau level, and the QP energies are close to the diagonal elements of the BdG Hamiltonian 



in the Landau levels representation, i.e. ||2J], [|T^]: 

E^^ = \j [Tiu, {n + 1/2 - nF)f + | A„,„ (J^) f (51) 
which is not an analytical function of Aq at the Fermi surface. This also explains the 



linear dependence of X{/S.Q,nF) on Aq, obtained by Norman et al. [jT3[ for small Aq at low 
temperatures. 

It is interesting to note that in our expansion the quadratic and the quartic terms for 
each Matzubara frequency v are proportional to qy and q^ respectively. Thus, the expansion 
parameter is actually x ~ 7- — A — go 

\TlUJc) Tip 
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In the high temperature hmit X > 1, where go ~ 2, it reduces to the temperature 
independent value x ~ A^/np used above. In the very low temperature limit ,X <^ 1, it 
diverges with (1/T)^, i.e. x ~ [Ao/vrfcsT'J^/ni? 

The breakdown of the small A expansion, resulting from this divergence at sufficiently 
low temperatures, seems to be related to the emergence of an opposite, high A expansion in 
the small parameter ^ ~ (vr/c^T/Ao)^, as obtained by Dukan and Tesanovic (see Eq.(|^)) 
in the low temperature limit. 

Application of the theory developed in the present paper to real experimental situations 
is not a straightforward matter; in addition to the influence of disorder in the vortex lattice 
and vortex fine fluctuations on the QP 'scattering' , discussed above, the 3D nature of the 
single electron band structure could also play an important role . The importance of the 
latter effect may be appreciated by noting that in contrast to the 2D model studied here, 
in a 3D electron system , e.g. with a spherical Fermi surface , Cooper-pairs in low Landau 
levels (i.e. for n,n ^ ) and with large longitudinal momenta ( i.e. near 
have the largest contribution to the SC condensation energy. This region is far away from 
the extremal orbit = 0, n = n^?, which dominates the dHvA oscillations . 

As a result, in addition to the QP near the extremal orbit, their counterparts with small 
cyclotron orbits (i.e. for n <^ np) and large longitudinal momenta kz , should also contribute 
significantly to the SC free energy in this case. The relatively strong sensitivty of QP with 
small cyclotron orbits to 'scattering' by the vortex lattice, as implied by the large damping 
parameter A found in Ref.. |13[ , may indicate that the QP 'scattering' effect in 3D systems 



is stronger than in the equivalent 2D systems. 

An effective parameter , np * < np , may be therefore introduced to take into account 
such an increase in the QP 'scattering' effect. 

Most of the SC materials in which clear dHvA oscillations were observed in the vortex 
state, such as V^Si, Nb^Sn, 1^X^2-82, and NbSe2, are essentially 3D systems with complex 
band structures and nonspherical Fermi- surfaces. One therefore expects characteristic values 
of n*p smaller than ni;' in these materials. Furthermore, the nonspherical Fermi surfaces , 
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combined with some unavoidable deviations from perfect crystaline order , should lead to 
some finite distribution around each dHvA frequency. 

A careful examination of the Fourier transformed spectra for these materials iiideed 
shows a significantly broad frequency distribution about each dHvA frequency F {=HnF for 
a spherical Fermi surface), with line width of the order AF ~ (100 — 200) T. This should be 
compared to the effective range of frequency modulation (AF).^^ = iJc2-f^mi)/4 {Hc2 — Hmv), 
associated with the expression within the square brackets in Eg. (^91) , which does not exceed 
15 T. Thus it is not surprising that the measured signal, which is the Fourier transform of 
this broad spectrum, does not exhibit a fine structure like the sign inversion predicted in 
our ideal 2D electron gas model. 

The organic superconductor k — {ET)^Cu {NCS)2 seems at first sight a good candi- 
date for testing the predictions of our theory, due to the quasi 2D nature of its electron 
band structure . Unfortunately, the transition from the normal to the SC state observed 
experimentally in this material is very broad p], extending far below the estimated value of 
ifmt,, which is found to be very close to the mean field value of Hc2 in this material. This is 
not surprising since the low dimensional nature of this compound and the low temperatures 
used in the dHvA experiments can lead to strong quantum fiucuations in the phase of the 
order parameter [^], [^], and so to the breakdown of the mean field approximation used 
in our theory. 

The relatively weak QP 'scattering' , predicted in the present paper, seems to be con- 
firmed, however, by the majority of the experiments performed so far: According to our 
theory it should lead to a significant deviation of the experimentally measured amplitude 
from the Maki-Stephen-Wasserman fitting formula (see Eq.(|l],H)) in the region where the 
leading SC effect exceeds the zeroth order (i.e. normal electron) term, i.e for H > Hin^. 
In this region the above qualitative analysis indicates that the damping of the dHvA os- 
cillations may be described by a parameter , A^/np * , smaller than the characteristic MS 
parameter A ~ Ag/y^np . Such a crossover from a relatively strong damping just below Hc2-i 
described well by the MS fitting formula, to a weaker damping at lower fields, was indeed 
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observed in almost all experiments carried out so far 

Furthermore, from the available experimental data two different characteristic slopes of 
the corresponding Dingle plot can be clearly distinguished . Our estimations show that the 
experimental crossover field Hcross from one slope to another is in a good agreement with 
the calculated inversion field Hin^- In particular, we have obtained for V^SilF = 1570T): 
Harass ~ 12.5r, H,nv ~ IS.ST, for YNi2B2{F = 51ir): Harass ~ 4.5T, H,nv ~ 4T, for 
NbSe2{F = 152T): Hcross ~ 5.6r, Hinv ~ 4T, and for NbsSniF = 581T): Hcross ~ 11.4T, 
Hinv ~ 13. 7T. 
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